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^ ■ Abstract 

' We consider the problem of network coding across multiple unicasts. We give, for wired and wireless 

O,' networks, efficient polynomial time algorithms for finding optimal network codes within the class of 

.^ , network codes restricted to XOR coding between pairs of flows. 

in 

■ I. Introduction 

In this paper we consider network coding across multiple unicasts, using the class of pairwise XOR 
codes introduced in [4] for wired networks. This class of codes includes "reverse carpooling" and two- 
^ ■ flow "star coding" for wireless networks [3]. We give efficient polynomial time algorithms for finding 

optimal network codes within this class on wired and wireless networks. 

> ' II. Problem setup 

Q> . A. Multiple unicasts problem on wired network 

T^lj- ■ K unicast sessions are transmitted over a network represented as a directed graph G = (A/", C) of 

O , N = \J\f\ nodes and M = \£\ links. Network coding is limited to XOR coding between pairs of uncoded 

Q . flows. Each session c = I, . . . ,K, demands a communication rate re- Each link (a, b) ^ C has a capacity 

^ I denoted Cab, which, if greater than maxcrc, is set to maXcTc. The total incoming capacity and total 

outgoing capacity of each node is upper bounded by C. A solution for a given set of demanded rates 

{vc} is an assignment of values to variables 
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We show how to find a solution for the problem with rates {rd if there exists a solution for the problem 
with slightly higher rates {(1 + 2e)rc} for any e > 0. 
We define a number of queues at each node i which can be interpreted as follows: 

• [/"': uncoded session c data previously at node v 

• p^''^ '^ : data from sessions c, c' coded at node j meant for both sinks 

• P^'^ •' : data from sessions c, c' coded at node j meant for sink dc 

• R'l'^ ■' : remedy for session c data that has been coded with c' data at node j 

Data can be transmitted on a link {a,b) E C from queue U^'" to Uj^^, from U^'" to C/^", from R'if-' to 
Rl'^ '^ , from pj^*^'*^ " to p^'^''^ '^ ^ or from Pa'^-' to -P^^'^''- A coding operation at a node a transforms / 
units from each of a pair of queues {U^"^', U^"" ) into / units in each of the queues i?^? "■ and i?^ ™. A 
decoding operation at a transforms / units from each of a pair of queues {Ra^-' ,Pa'^ ■') into / units in 
C/a"'. A branching operation at a transforms / units from a queue p^'^''^ " into / units in each of the 
queues Pa*^ -^ and Pa '^■' . 

Then the solution variables can be interpreted as follows: 

• v"^: average flow rate from U™ to C/™ 

• v"^""- average flow rate from U™ to U^"" 

• '^ab'^ '■ average flow rate from Pa to p^'^''^ " 

• vr^^ •' : average flow rate from Pa'^ •' to P^^'^ •' 

• p'^fjj -' : average flow rate from R'if •' to R'^^^ -^ 

• Ja''^ '" ■ average rate of coding transformation from (f/™, U^^'^') to (-R^? ", i?^ ™) 

• r]a'^-': average rate of decoding from {R'^a\Pa^'') to [/«"' 

• a a '■ average rate of branching transformation from Pa''^ to {Pa'^ ■' ,Pa'^-') 

B. Modified problem 

For a given problem instance, we consider a modified problem any solution of which is equivalent to 
a solution of the original problem and vice versa. The modified problem reverses the direction of the 
poison flows and imposes some additional constraints. Specifically, it is defined as follows: 

• Each source node s^ has a source queue U'^, an overflow queue tJ'^, and a virtual source link of 
capacity C from V^ to queue 11^'^'=. 

• Each node a has a virtual coding link, a virtual decoding link and a virtual branching link, each of 
capacity (7/2. 

• Each real and virtual link e is associated with a set Ve of pairs {O, V) such that data units from 
each queue in the set O are transformed into an equivalent number of units in each queue in the 
set V via e: 

- if e is the virtual source link for session c, Ve = {U^, U^^") 

- if e is a real link (a, b) G C, 

Ve = {(f/r , un, {K\ un, {pt''^\pt''^^), (p^'^ ^r'^), (i?f ^ c^h 

- if e is the virtual coding link at node a, Ve = {{{U^" , U^"'), {R^?''', R^""")) ■ Cy^c';a^v, v'} 

- if e is the virtual decoding link at node a, Ve = \ [ R'a '' , [Pa^ "' , Ua'' ) ) : c^ c';a ^ j\ 



- if e is the virtual branching Unk at a, Ve = \ i {Pa^-' ,Pa ^■'),Pa ) : c / c'; a / j > 
• Data is removed from queues [7™, P'^ '^'' , and p^'^''^ '\ 

C. Wireless case 

We also consider the multiple unicasts problem and corresponding modified problem in a wireless 
setting, using the following wireless network model. We model wireless transmissions by generalized 
links, denoted by (a, Z), where a is the originating node and Z is the set of destination nodes. The 
network connectivity and link transmission rates depend on the transmitted signal and interference powers 
according to some underlying physical layer model. For example, the transmission rate per unit bandwidth 
fiij from node i to node j, with other nodes n ^ J\f transmitting independent information simultaneously, 
may be given by the Shannon formula [2] 

/..,(£, S)= log fl+ ^'^'^ 



where Pi is the power transmitted by node /, Sij is the channel gain from node I to node j and A'^o is 
additive white Gaussian noise power over the signaling bandwidth. For simplicity, we consider a finite 
set C of links and a finite set U of sets of simultaneously achievable link rates. We denote by C(^a,z),u 
the capacity of link (a, Z) in set u ^U. K solution to the multiple unicasts problem consists of a convex 
combination of sets in lA, which gives a set of average link capacities achievable by timesharing, and a 
network code that operates over the network with these average link capacities. 

In the wireless case, each node has a virtual source link, coding link and decoding link but no virtual 
branching link; a branching operation occurs over a real wireless link. For virtual source, coding and 
decoding links e, the set Ve is defined exactly as in the wired case. For real wireless links (a, Z), 

na,z) = {{u^\ un, {K\ un, {pt''^\pt'^'), {pf'\pf'), {Ra'^^K'n, 

((t/-, ijf ^■), (C/r, 4"')), {{Pf\ P'^,'''),Pt'^^) : c + d- b, b' G Z} 

III. Back-pressure approximation algorithm 

Queues U™ , V^, R'j^'" , Pa'^ '' are associated with session c. We consider each joint poison queue pj'^''^ " 
as a pair of queues, one associated with each session c and c', such that equal amounts are always added 
or removed from each of the pair of queues. 

The overflow queue for each session c has a potential aje"'"^^" that is a function of its length /, where 

"^ = 24k' ^'^ 

and B is a. constant to be determined. We divide each other queue into subqueues, one for each link for 
which it is an origin or destination. Each subqueue Q has a potential 't>c{Q){^Q) = e^'^'Q''" that depends 
on its length Iq and its session c{Q). For notational simplicity, we will abbreviate subscripts c{Q) as 
subscripts Q. 
Flow entering or leaving subqueues associated with session c is partitioned into packets of size pc = 
(1 + e)rc. Besides its true length Iq, each subqueue Q has an approximate length Iq that is an integer 
multiple of the packet size. The approximate length of a subqueue is updated only when its true length 



has changed by at least one packet since the last update of its approximate length, as follows: Iq is set 
to kpQ where k = -^^— if Q is an origin subqueue or k = \ -^^ if Q is a destination subqueue. 
Between updates, Iq and Iq satisfy Iq — 3pQ < Iq < Iq for an origin subqueue, ox Iq < Iq < Iq + . 
for a destination subqueue. 
We denote by Pe the subset of Ve consisting of pairs (O, V) € Ve satisfying 

mill /o > 

maxTg < Brg + ln((L + l)p)/aQ + 



where 



maxc Tc 

P= — 

mm^ Tn 



A. Wired case 
In each round t, the algorithm carries out the following: 

1) Add (1 + e)rc units to the overflow queue tj'^ of each session c, then transfer as much as possible 
to V^ subject to a maximum length constraint of Brc for V^ 

2) For each real and virtual link e, flow is pushed for zero or more origin-destination pairs (O, V) G Ve 
such that the total amount pushed is at most the link capacity Cg. Specifically, initialize C to Cg 
and repeat 

• Choose the pair (O, V) G Pf. that maximizes 

Let 



C' = min(c, ^mm^^(pQ-|<^Q|)) 



where 5q is the change in Iq since the last update of Iq. Subtract C units from C, subtract 
C units from Iq for each Q € O and add units C" to Iq for each Q € I?. For Q € (O U V), 
if C' = pQ - \5q\, update Iq. 
. If C = then end. 

3) Zero out all subqueues Uf^ , P^ '^■' and P^'^'^ ". 

4) For each queue that has at least one subqueue whose actual length has changed during the round, 
reallocate data units to equalize the actual lengths of all its subqueues. If the actual length of any 
subqueue has changed by at least one packet since the last update of its approximate length, update 
its approximate length. 

When the amount of flow remaining in the network queues is an e-fraction of the total amount that has 
entered the network, the flow values for each link are averaged over all rounds to give the solution. 



B. Wireless case 

The algorithm is the same as for the wired case, except for Phase 2, which is as follows. For each virtual 
link e, flow is pushed for zero or more origin-destination pairs (O, V) € Ve exactly as in the wired case. 
For each real link e = {a,Z), flow is pushed for zero or more origin-destination pairs {0,V) G Ve 
such that the total amount pushed is at most the average link capacity XuCe,u for some A„ satisfying 

Ylueu ^u < 1- 
Specifically, initialize T to 1 and repeat 

• For each real link e = (a, Z), let 

We = max_ W(o,v) 
{Oe,'De) = arg max w/o,v) 
where wiq^j)\ is defined in ^. Choose the set u gU which maximizes 

Vu = ^ WeCe,u- (3) 

e 

Let 

rrf ■ frr ■ ■ PQ-\h\ 
I = mm I , mm mm — — 

V e:Ce,„>OQe{C'.U©e) Ce,u 

where 6q is the change in Iq since the last update of Iq. Subtract T' units from T, and for each e, 
subtract T'Ce,u units from Iq for each Q ^ Oe and add the same amount to Iq for each Q € Ve, 
updating Iq if T'Ce,u = PQ - \Sq\- 
. If T = then end. 

IV. Potential analysis 

Our analysis is partly based on the approach in [1]. We lower bound the decrease in potential over a 
round t. We denote by Q{t) the actual length of a subqueue Q at the end of Phase 1 of round t. From [1], 
the increase in potential during Phase 1 is upper bounded by 



^Pe<P'eiU%t)). (4) 



A. Flow solution-based algorithm 

We lower bound the decrease in potential during Phases 2 and 3 by comparison with the potential 
decrease resulting from pushing flow based on a flow solution for rates /c = (1 + 2e)rc. This algorithm 
differs from the back pressure algorithm only in the specific portion of Phase 2 determining the amount 
of flow pushed for each pair (O, V) G Ve of each link e. 

Consider a flow solution for rates fc- Partition the flow for each session c into elementary flows T^ such 
that all data in an elementary flow undergoes the same routing and coding operations. Each elementary 
flow T^ has a size denoted /^, and consists of a set of links with associated origin and destination 
subqueues, comprising a primary path from Sc to dc and a remedy path associated with each coding 
node. Note that each flow T^ starts from queue W^, which consists of a single subqueue, and that 



fc = X^^ fn- Let L be the length of the longest primary path and F the maximum number of links in an 
elementary flow. A subqueue in T^^ is considered upstream or downstream of another according to the 
direction of flow in the modified problem defined in Section III-BI 

Phase 2 of the flow solution-based algorithm consists of a preprocessing part and a flow pushing part. 

1) Preprocessing procedure: 

• Initialization: 

Remove from each ^^ any portions of each branch that are downstream of a subqueue Q in T'f^ 
for which Iq < 3pc- Note that all subqueues of each queue have been equalized in Phase 4 of 
the previous round. Let Q be the set of flows J^^ containing some subqueue of length at least 
Brc + ln{{L + l)p)/ac- For flows T^ G Q, let Qn be the furthest downstream origin subqueue in JF^ 
of length at least Br^ + ln((L + l)p)/ac', for flows T^ ^ Q, let QJi be the longest origin subqueue 
in JF^ (if there is a tie, choose the furthest downstream). 

Initialize H as the set Hq of all flows JF*"^, initialize Hi and H2 as empty sets, and for each c, n, 
initialize Q*l^ as Q^- 

• Phase A: Repeat 

- Choose some flow T^ (^ Q nH and remove from H all flows J^^, such that T^, shares any 
portion of a common joint poison path segment, or the coding and branching links at either 
end, with T^. Remove J^^ from Tt and add it to Hi. 

- If there is no such flow remaining, end Phase A. 

• Phase B: Associate with each flow JF^ G 7Y a weight w^, initially set to (/)^([/'^(t)). For each flow 
J^^, G H, set T^, to be the set of subqueues Q'^ such that T^, is the shared flow corresponding to 
Q^, and either 

- Q^ is a remedy or individual poison subqueue whose corresponding branching link is in JF^ 
and whose corresponding coding link is in T^,, or 

- Q^ is a joint poison subqueue whose corresponding branching link is in ^^,. 
Repeat 

- Choose some flow JFi;, G H such that 



E 



<:< V (<A^(Q;i(t))-o, (5) 



* remove T^, from H, 

* remove any subqueues in .F^, from T^ V c, n, 

* set the weight w'^ of each flow JF^ G 2^, to (/)^(Q^(t))- 

- If there is no such flow, for each 2^, and each Q'!^ G T^/ set Q*"^ to C/^, and end Phase B. 
Phase C: Repeat 

- Choose some flow T^ ^ H such that the longest individual poison subqueue along one of its 
poison branches, Q, is longer than Q*^, and the entire joint portion of that branch together with 
the branching link are not in the corresponding shared flow. Set Q*^ to Q and remove from T^ 
all but the portion downstream of Q. 

- If there is no such flow remaining, end Phase C. 



Observe that: 

• At most L flows are removed by each flow in Q, and 

0;(Sre + ln((L+l)p)K) > (L+l)max(/.^(5re) 

c 

> (L + l)max(/)^(C/'=(t) (6) 

c 

• fn = fn' for aU ^^ and JF^', that share a common joint poison path segment. 

• At the end of the preprocessing procedure, 

2 E ^'ciQnit))fn > Yl ^'ciQnit))fn- (7) 

To see this, note that at the end of Phase B there is a one-to-one correspondence between flows 
Tn^H for which Q*^ / Q'^ and elements Q^ in the sets Z^/ of flows J^^' ^ ^- For ^^^h J?^^, G H, 
by (H, 

<''> E {<t^'c{Qnit))-0- 

Multiplying by /^ and summing over all T"^, G Ti, and noting that at the end of the preprocessing 
procedure 

(t^'Mnit)) > < (8) 

gives (0. 

• Phase B and C maintain the invariant 

E <fn > E '^c(C^^(t))/n (9) 

where H2 is the value of H at the start of Phase B. The invariant holds since both sides are equal 
at the start of Phase B, and the left-hand side is monotonically non-decreasing. From ((SJl, © and 
^, at the end of the preprocessing procedure we have 

E ^'ciQnit))fn> E ^'ciU%t))f^. (10) 

2) Flow pushing procedure: For each flow JF^ G TYi U H, let T*'^ be the portion of T^ downstream of 
Q*"^. Partition its links into a set Q^'^ of subsets. Each of these subsets S G Q^'^ may be 

• a path from Q"!^, if it is a poison subqueue, to its associated coding node a{S) 

• a path from Q"!^, if it is a remedy subqueue, up to and including its associated decoding link at 
node w{S), and the associated poison path from w{S) to its associated coding node a{S) via the 
branching hnk at node b{S), 

• a path associated with uncoded flow ending in sink node dc, or 

• a path associated with uncoded flow ending in a coding link at a node a{S), together with the 
associated remedy path up to and including the decoding link at a node w{S), and the associated 
poison path from w to a via the branching link at a node b. 

Note that each subset starts either at Q^^ or at an uncoded subqueue. 



The flow solution-based algorithm pushes flow as follows. First, for each pair {Sc,Sc') € Q*'^ x Q*^? 
that shares any portion of a common joint poison path segment, note that /^ = /^, and that Phase B of 
the preprocessing procedure ensures that both subsets or neither contains the coding link; in the latter 
case both or neither contains the branching link. Thus, the following are the only two cases: 

• Case 1: One of the subsets, say Sc, contains both the coding link at a{Sc) and the branching link at 
b{Sc), while Sc' contains the coding link but not the branching link. Phase C of the preprocessing 
procedure ensures that all session c individual poison subqueues along the joint poison path segment 
are shorter than Q'^^. Push /^ units through ScUSc', pushing session c individual poison units through 
the joint poison path segment. 

• Case 2: Both subsets contain the same portion of the joint poison path segment. Push /^ units 
through Sc U Sc' ■ 

Next, for each subset S of some J^*^ that does not have any coded segments in common with T*? for 
all c' 7^ c, n', we have the following cases: 

• Case 1: Q^ is an individual poison subqueue in S. Phase B of the preprocessing procedure ensures 
that Q^^ is the longest individual poison subqueue in J^*?' . Push /^ individual poison units through 
S. 

• Case 2: Q"^ is a joint poison subqueue in S. Push /^ joint poison units along the path in S 
downstream of Q**^. 

• Case 3: Q'^ is a remedy subqueue in S. Phase C of the preprocessing procedure ensures that Q"^ is 
longer than all session c individual poison subqueues along the primary path of S. Push /^ remedy 
units along S through the decoding link at w{S), and /^ individual poison units along the primary 
path of S. 

• Case 4: Q*^ is an uncoded subqueue or is not in S. Push /^ uncoded units along the primary path 
of S starting from its longest session c uncoded subqueue. 

Note that flow is pushed only from origin subqueues Q for which 

Iq > 3PQ ^ ^Q > 
and only to destination subqueues Q for which 

Iq < Brg + ln((L + l)p)/aQ => Iq < Brg + ln((L + l)p)/aQ + Spg. 

B. Potential decrease in Phases 2 and 3 

The decrease in potential from pushing / units across a link from a set O of origin to a set T> of 
destination subqueues is at least 

where Iq denotes the initial length of each subqueue Q. 



1) Flow solution-based algorithm: We denote by 0(^*^) and V{T*'^) the sets of origin and destination 
subqueues of a flow JF*"^, and by Oc,e and Pc,e the sets of session c origin and destination subqueues of 
a link e. For each subqueue Q, denote by fg the total flow out of Q (if Q is an origin subqueue) or into 
Q (if Q is a destination subqueue) in the flow pushing procedure of the flow-solution based algorithm. 
The potential drop over Phases 2 and 3 in the flow solution-based algorithm is at least 

El E ifQ^'cih)-fQ€ilQ + fQ))- E fQ{<P'cih)+fQ<Pc{lQ + fQ)) 
^ El E fQ{<t>'cilQ)-Mc{lQ + fc))- E fQ{<P'cilQ) + fc<PcilQ + fc)) 

> El E fnWc{lQ)-M'c{lQ+fc))- E fnmQ)+fc€{lQ + fc)) 

Since all subqueues of each queue have been equalized in Phase 4 of the previous round, and since each 
flow T*^^ has one origin subqueue of length Q*j^{t), at most L + 1 destination subqueues each of length 
less than 3pc> and a total of at most F links, this potential drop is lower-bounded by 

E/n i^'ciQnit)) -{L+l) cP'.iSp,)) -Y^Frj,<t>c{Qn{t) + fc) - Y. ^F fnfcCt>'^{Qn{t) + fc) 

From Q, we have 

3FW^{l + f,) < 3F(l + 2e)r,0:f(/ + (l + 2e)r,) 
= 3F(1 + 2e)rca^e"^'+°<=(i+2^)^-= 
= 3F(1 + 2e)r,ace"=(i+2eK^/^(^) 

= ^ e 24F 0^(/) 

< e</'c(0/4. (11) 

This yields, using 0, the following lower bound on the potential drop: 

YfU^'ciQnit))-iL + l)<P'ci3pc))-Yyn<^<P'ciQn{t))/4- E fne<Pc{Qnit))/^ 

> E /n ((l - i) '/'c(Qr(t)) - (^ + 1) -^cCSPc) 

> 5^(1 + 26)r, ( (l - I) 0^(C/^(t)) - (L + 1) 0;(3pe) 

c 

> E^^((l + |-^')'^c(^'(*))-(^+l)(l + 2e)'^c(3Pc)) (12) 
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2) Back pressure algorithm: If the packet size Pc were infinitesimally small rather than (1 + e)rc, 
then Iq = Iq and the procedure in Phase 2 of the back pressure algorithm would give, for each link 
e, the maximum possible potential decrease from pushing flow for zero or more origin-destination pairs 
(0,P) G Pe such that the total amount pushed is at most the link capacity Cg. Since pc = 9(rc) and 

0^(/ + e(r,)) - 0^(/) < Q{r,)^:{l + e(r,)), 

the back pressure algorithm achieves, for each link e, a potential decrease of at least that achieved by 
the flow solution-based algorithm minus an error term 

QG(au©e) 

where fq is the amount of flow added or removed from Q in Phase 2 of the flow solution-based algorithm. 
Since /q = G(rQ), decreasing each Qc by some constant factor is sufficient to ensure that (fT2t applies 
to the back pressure algorithm. 

C. Overall potential change and number of rounds 

The potential does not increase during Phase 4. Thus, from Q and M2\ . the overall potential decrease 
during the round is lower bounded by 

E ^^ ((I - ^') ^'ciU%t)) - (L + 1) (1 + 26)<^^(3p,)) 

c 

= E ^= (i - '') '^c(^'W) - (i + 1) (1 + 2^)^^"^ 

c 

If t/'^(t) = -Brc for some c, then the decrease in potential is at least 



which is non-negative if 



r,a, [^ - e') e^^^^^ - K {L + 1) (1 + 2e)^^e— 

If U^{t) < Bvc for all c, the overflow queues are empty and, since there are Q{NMK) session-c 
subqueues each of potential less than 

(t)c{Brc + ln((L + l)p)/ac + Ap^) = 6 [ipe^^ , 

the overall potential in the system at the end of the round is at most 

e (NMK'^Lpe^^ 
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By induction, this is also an upper bound on the total potential at the end of every round. The length of 
the overflow queue for session c is thus never more than 

f NMK^Lpe^\ _ fNMK^Lp' 



/NMK^Lp\ 



and the total units for session c is at most 

0(^^^^^ + iVMK(i?., + ln(Lp)K + 4p,)) = 0(^(KL, + ln(^) 



Thus, at most 






rounds of input flow for each session remain in the network at any time. For 

amount remaining in the network is at most a fraction e of the total amount that has entered the network 
up to round t. 

D. Number of operations 

In each round, at each node, pushing flow across links results in a total decrease of at most 0{C) 
in the actual lengths of origin subqueues and a total increase of at most 0{C) in the actual lengths of 
destination subqueues. The total change in subqueue lengths from rebalancing is not more than the total 
change resulting from pushing flow across links. Thus, at most 0{ND) approximate subqueue lengths 
in the network are updated in each round, where D = maxc C /re- Assuming the subqueue differences 
computed by the algorithm can be stored, only those differences involving subqueues whose approximate 
lengths have changed are recomputed. 

For the coding link at each node i, for each pair of sessions {c^d}, the values of wi^q^x>)^ defined in 
©, for pairs {0,V) G Pe of the form ((C/f^, C/f"'), {R^/^R^^"')) are stored in a sorted list of length 
0{N'^), and the list indexes are ranked according to the maximum value in each list. If U^"" or Uf" 
changes, A^ differences are updated in each of K lists, requiring 0{NK log N) operations, and the rank 
of these K lists are updated using 0{KlogK) operations. If i?^f * or R'^'^^ changes, N differences are 
updated in one of the lists and fewer operations are required. 

For the decoding link e at each node i, the values of w^qjj^ for pairs (O, V) E Pe are stored in a 
sorted list of length 0{NK'^). A change to any of the approximate subqueue lengths requires updating 
of at most K differences, for which 0{K\og{NK)) operations suffices. 

1) Wired case: For each real and virtual branching link, the values of wtQm for pairs (O, V) G P^. 
are similarly stored in a sorted list of length 0{NK'^). A change to any of the approximate subqueue 
lengths requires updating of 0(1) differences, for which 0{\og{NK)) operations suffices. 

Thus, each round has complexity 0{N'^KD\og{NK)), and the algorithm has complexity 



12 



2) Wireless case: We assume that for each wireless Unk (a,Z), \Z\ = 0(1). For a wireless link 
e = {a,Z), for each pair of sessions {c,c'}, the values of ^(o©) for pairs {0,T>) G Ve of the form 
(([/„'', Ra"^), {U^""-, Rb'"^)) are stored in a sorted list of length 0{N'^), and the values of W(^a.v) for all 
other pairs {0,T>) € 'Pe are stored in a sorted list of length 0{NK'^). The values of y^, defined in ^, 
for each set n G ^ are also stored. 

Of all the subqueues of {a,Z), a change in the approximate length of some subqueue U^'" or C/^*^ 
requires the most operations: 0{N) differences are updated in each of 0{K) length-0(A^^) lists, 0(1) 
differences are updated in the length-O(A^-fr^) list and 0{K) of the 0{K'^) list indexes are repositioned, 
requiring a total of 0{NK\ogN) + 0{\og{NK)) + 0{K log K) = 0{NK\og{NK)) operations. A 
further 0(|^|) operations suffices to update the values of yu and find the maximum among them. 

Thus, each round has complexity 0{ND{NK\og{NK) + \U\), and the algorithm has complexity 

P ^NH,KFDiNK^o^(NK) + m (^„^^ ^ ,„ |^^j^ j 
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